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Abstract - We evoke the idea of representation of the chaotic attractor by the set of unstable pe- 
riodic orbits and disclose a novel noise-induced ordering phenomenon. For long unstable periodic 
orbits forming the strange attractor the weights (or natural measure) is generally highly inhomo- 
geneous over the set, either diminishing or enhancing the contribution of these orbits into system 
dynamics. We show analytically and numerically a weak noise to reduce this inhomogeneity and, 
additionally to obvious perturbing impact, make a regularizing influence on the chaotic dynamics. 
This universal effect is rooted into the nature of deterministic chaos. 
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Introduction. — During few recent decades, noise 
was found to have ability not only to make an obvious 
disordering impact, but also to play a constructive role, 
increasing degree of order in system dynamics P~HI] ■ Most 
remarkable advances are the phenomena of stochastic [T] 
and coherence resonances [2J; the suppression of deter- 
ministic chaos in dynamics of dissipative systems by noise 
(thoroughly discussed, e.g., in [3]); synchronization by 
common noise [I] ; Anderson localization [3] ; the analog of 
the latter phenomenon for dissipative systems [5J; etc. In 
this Letter we disclose a novel phenomenon, where weak 
noise reduces disorder in chaotic dynamics. The effect 
is rooted in inhomogeneity of the distribution of natural 
measure over the strange attractor and impact of noise on 
this distribution. 

First, we demonstrate this inhomogeneity and reason 
for the consequences expected from it. Then we construct 
an analytical theory for the noise impact. Finally, the 
predictions of the analytical theory arc underpinned with 
results of simulation for paradigmatic chaotic system, the 
Lorcnz one [9], where one can see the disorder reduction 
effect with divers quantifiers of regularity, and for the map- 
ping which is relevant for one of the classical routes of the 
transition from regular dynamics to chaotic one [8]. 

In our treatment we relay on the idea of representation 
of chaotic attractor by the set of UPOs embedded into it ( 
[TO j lll j . etc.). On this way, for instance, one can evaluate 
the average A for the chaotic regime as A = J A(x)n(x)dx, 
where fj, is the natural measure, and integration is per- 
formed over the whole attractor. It was shown in |10j 



(see also the text book [TT|) that the distribution of natu- 
ral measure over attractor can be recovered from proper- 
ties of the UPOs. For two-dimensional maps (and there- 
fore three-dimensional dynamic systems where one can use 
Poincare section to construct such a map) the natural mea- 
sure of an UPO is fi oc A -1 where A is the largest multi- 
plier of perturbations of this UPO. This result are treated 
to be valid in higher dimensions as long as we have only 
one positive Lyapunov exponent. The approach was em- 
ployed already in [T2] for evaluation of fractal characteris- 
tics of the Lorenz attractor from the UPOs data weighted 
with measure /i oc A -1 . 

Inhomogeneity of natural measure distribution. 

— We consider the Lorenz system as an example, 



dx dy 
— = <j(y — x), — 
dt yy h dt 



dz 

y-xz, — = xy - bz. (I) 



Recently [13], we reported the results of calculation of 
all UPOs with length N < 20 (by "lenght" of UPO we 
mean the number of loops). In Fig Q] one can see proba- 
bility density function (PDF) of (z) calculated over UPOs 
of length N — 20, the distribution weighted with A -1 
and PDF of (z) calculated over 20-loop segments of the 
chaotic trajectory. The unweighted and weighted distri- 
butions are remarkably dissimilar and the latter matches 
the chaotic averaging much better. The weighted distri- 
bution possesses a relatively big tail for large (z), where 
there is a quite few number of UPOs (unweighted PDF 
nearly riches zero next to (z) = 24). The origin of this dif- 
ference is a huge dispersion of multipliers for long UPOs: 
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Fig. 1: (Color online) Probability density function of (z) cal- 
culated for the Lorenz system with classical set of parameters 
O = 10, b = 8/3, r = 28) over UPOs of length N = 20 with- 
out weight (doted green line), with weight A -1 (solid red line), 
which corresponds to UPO's measure within chaotic attractor 
in noise- free case, and with weight VA/A (dashed blue line), 
which corresponds to the maximal correction of the weight in- 
duced by a weak noise, and 20-loop segments of chaotic orbit 
(solid black line). 

for N = 20 the smallest multiplier A = 62 200 is by factor 
72.8 smaller than the largest one A = 4 527 160. Notewor- 
thy, the weighted distribution is significantly broader than 
the unweighted one. This feature implies a large contri- 
bution of relatively small number of UPOs with extreme 
deviation of properties from the typical ones (e.g., in Fig.Q] 
these UPOs have untypically large values of (z)). In terms 
of temporal evolution, the system dynamics is strongly in- 
homogcncous: the epoches of "typical" behaviour intermit 
with the epoches spent near several least unstable "untyp- 
ical" periodic orbits. 

Such a huge inhomogcncity of multipliers evokes the 
questions: (i) How this inhomogeneity is effected by weak 
noise? (ii) What are the consequences for the system dy- 
namics? Indeed, the set of UPOs is dense and even in- 
finitely small perturbation throws the system from one 
UPO onto another. Hence, one can expect that a very 
weak noise could first influence the rules of walking over 
this set and only then affect properties of UPOs. Pre- 
sumably, the noise should play a "smoothing" role, re- 
ducing inhomogeneity of multipliers. As a result of such 
smoothing, one should observe, in particular, shrinking of 
the distribution of (z) calculated over finite segments of 
the chaotic trajectory because unweighted distribution is 
more narrow than the weighted one. One could also ex- 
pect other kinds of rcgularization of the chaotic dynamics 
by weak noise. 

Effect of noise on natural measure distribu- 
tion. — A significant advance in analytical description 
of noise action on systems with fine structure was made 
for one-dimensional maps [HI [15]. In [14] the formal- 
ism of Feynman diagrams was employed for treatment 
of noisy dynamics as a perturbation to deterministic dy- 
namics over the set of periodic orbits. In [15] the escape 




Fig. 2: (Color online) Cell on the Poincare surface. In the 
vicinity of fixed point A of map F fe , curvilinear parallelogram 
B'C'D'E' is the image of BCDE. 

rates for unimodal maps were evaluated. However, in 
higher dimensions — which are practically inevitable when 
the map is a Poincare map of chaotic continuous-time 
dynamic system — analytical advance becomes restricted; 
studies require ad hoc approaches and/or more numer- 
ics [7J. Moreover, for most problems addressed in liter- 
ature, low-length periodic orbits provide enough data for 
an accurate description [141115] , whereas in our case inho- 
mogeneity of high-length orbits is essential. 

For evaluation of natural measure, we first briefly recall 
its calculation for the no-noise case (interested readers can 
find a detailed rigorous consideration in Sec. VI of [TO]). 
We consider certain Poincare section and two-dimensional 
map on it, which is induced by three-dimensional phase 
flow. We are interested in calculation of the natural mea- 
sure distribution over attractor on the Poincare surface. 
Let us consider k iterations of the map, x„ +/ t = F fc (x n ), 
where k is large, and a small cell bounded by stable and 
unstable manifolds, which feature the directions of com- 
paction and stretching of the phase volume, respectively, 
as sketched in Fig. [2] In the vicinity of fixed point A of 
map F fe , curvilinear parallelogram B'C'D'E' is formed by 
states belonging the cell before k iterations of the map, 
and BCDE is the pre-image of this parallelogram. By 
definition of the ergodic (mixing) attractor, for any two 
subsets S a and S& in the phase space of the system, we 
have 

lim [i{S a r)F k (S b )} = n{S a ) fi(S b ). 

k— too ' ' 

If we take S a and Sb to be the cell we consider, then 
/i(cell) = //(BCDE)//i(cell) = /i(B'C'D'E')/^(ccll). Fur- 
thermore, the measure varies smoothly along the unsta- 
ble direction [16] (this is quite general in nature, that 
the stretching/expansion is a smoothing process), while 
it typically has a complex inhomogeneous structure in the 
stable direction, which is transversal to the chaotic attrac- 
tor. Therefore, one can easily compare the measure of the 
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cell and parallelogram BCDE: their ratio is the ratio of 
their extension along the unstable direction, because these 
two sets have one and the same structure and extension 
along the "complex" inhomogeneous stable direction. We 
find //(cell) = ,u(BCDE)//j(celi) = CD/CD' = l/A pk (A), 
where Apt (A) is the largest eigen value of fixed point A 
of map F k . Recall, that we have found the measure as- 
sociated with a single fixed point A; when there is more 
fixed points in the cell one has to sum contributions of all 
of them, 



/x(cell) = lim 



E 



i 



fixed points in cell 



A F *(Aj) 



Henceforth, the subscript j indicates the number of a fixed 
point or the corresponding UPO. This result was derived 
in [TU]. Recall, the presented evaluation implies that at- 
tractor is hyperbolic, i.e., the stable and unstable man- 
ifolds are not mutually tangent anywhere. Although the 
Lorenz attractor is non- hyperbolic, it is generally accepted 
that statistically significant measure distribution is deter- 
mined by the same law 1/A. 

This calculation of the natural measure sheds light on 
the fact that we have to track only the coordinate along 
the unstable manifolds, say 77, and evaluate the fraction 
of states in the e-vicinity of some UPO, for which 77 stays 
within the e-vicinity. In no-noise case, the fraction of such 
states is simply 1/A. Now let us evaluate such a fraction in 
the presence of a weak Gaussian 5-correlated noise a£(t): 
= 0, (£(*)£(*')) = 26(t-f), a is the noise amplitude. 

Our map is induced by a continuous time evolution and 
in the presence of time-dependent signal, noise, one has to 
deal with a continuous-time evolution. Placing the origin 
of the ry-axis at fixed point Aj , one can find 



+ a /[x,(T)R(t +T)e A ^- T >d7 



(2) 



where Tj is the period of UPO Xj(i) associated with Aj, 
/(x) is the susceptibility to noise, Aj - Tj~ kiAj is the 
leading Lyapunov exponent of the j'-th UPO, R is the 
Gaussian random value of unit variance, (R 2 ) — 1, and 



2 _ 2 
a = a 




2a 2 / 2 [x,(r)] e^^^dr = ^ 



A, 



Here Fj is the characteristic value of /[xj(r)]; F 2 = 

2\ 3 (l-Aj 2 )- 1 / Tj / 2 [x 3 -(t)] exp(-2A i T)dr. According to 
Eq. J2]), the probability of r)(to+Tj) to be at the ej-vicinity 



of fixed point Aj is 

P(|»?(*o+2j)| <ej\v(to)) 



-'ii('o) 



2tt 



dR 



~ E j _A ? T ?( t o) 



<I> 



Ej - AjTyfa) 
ay/2 



aV2 



where $(...) is the error function. Finally, the measure of 
the £j -vicinity of fixed point Aj, say fl Sj (Aj), is 



(3) 



Aj [FjaM 2 



The last approximate expression corresponds to the limit 
of large Aj, which is relevant for long UPOs we consider. 
For vanishing noise, a = 0, the error function turns 1, and 
we have the conventional result of [10] , 

One can see in Eq. (J3j> , that for larger Aj (hence, larger 
Aj) the argument of $ is larger. This decreases inequiva- 
lence between UPOs with different weight. Maximal de- 
crease of inequivalence is achieved when the argument 
of $ is small, (£j/Fj)y/\/2 < o; at this case one finds 
/i oc ^/Aj/Aj. In particular, for the 20-loop UPOs in the 
Lorenz system the maximal ratio of weights decreases by 
15%. In Fig. Q] the distribution of (z) with weig ht v^/A 
is slightly shrunk compared to the one with A -1 . 

Notice, the vicinity extension e is not precisely known 
value in our analytical theory The properties of neighbor- 
ing UPOs could be similar and, additionally, the measure 
smoothly varies along the unstable direction on the attrac- 
tor. Hence, e could be not only the characteristic distance 
between nearest UPOs of length N but rather a consider- 
ably larger number featuring the size of clusters of UPOs. 

Effect of noise on disorder quantifiers. Now we 

calculate the variance of (z) for UPOs of length N in the 
presence of weak noise a£(t): 



var«z)j) = J2 Pi [( z fj - < z > 2 ] + K( 



(z) 



j 



(4) 

Here we sum the variance related to the distribution over 
the set of UPOs of the noise-free system and the noise- 
induced distortion of these UPOs themselves, which is 
nearly statistically independent from the former and de- 
scribed by term Ka 2 . The value of parameter K is approx- 
imately inferred from numerical calculation of the variance 
of (z) for the segments of the chaotic trajectory (rough 
analytical estimation yields the value of the same order of 
magnitude). 

In Fig. [31 we plot the variance of (z) for UPOs of length 
N = 20 for the Lorenz system ([T]) with the classical set 
of parameters. We approximately assume Sj/Fj to be the 
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Fig. 3: (Color online) Theoretical dependence of the variance 
of (z) for UPOs of length N — 20 on noise amplitude a for 
the Lorenz system with classical set of parameters and various 
values of parameter e which measures characteristic vicinity of 
UPOs. For the sake of demonstration we plot the first term 
of the expression for the variance [Eq.@] (solid lines) which 
features purely noise-induced redistribution of measure over 
UPOs of the noise-free system; the dashed lines with symbols 
represent the total variance with account for the noise-induced 
dispersion of each single UPO. 

same for all UPOs (notice, Fj is of order of 1). One can 
note the always-present ordering role of noise which de- 
creases inequality of weights of UPOs: solid lines show 
decrease of the variance by 2.5%. However, this order- 
ing effect can be overwhelmed by the dispersive action of 
noise on the orbits when e is not small enough. For the 
classical set of parameters the effect is not enough well 
pronounced and practically unobscrvablc with dispersion 
of (z) for chaotic segments. Still, we report the results for 
this case, where the effect is small, for to emphasize their 
universality: the effect is always present in some form. For 
a stronger inhomogeneity of the attractor one can expect 
a more pronounced effect. 

We performed simulations for the Lorenz system ([T]) 
with additive noise a£(i) put in dz/dt. With the classical 
values of parameters a = 10 and b = 8/3, the inhomogene- 
ity of the chaotic set is strongest near the threshold where 
it becomes attracting, r ps 24.06 (e.g., see [17j[T8]). At 
r ps 24.74 the chaotic attractor becomes the only attrac- 
tor in the phase space (between the stability threshold and 
this point, it coexists with the stable fixed points which 
significantly influence the noise-perturbed dynamics of the 
system). In Fig.HJa), the dependence of var(z) for 20- loop 
segments of the chaotic trajectory on the noise amplitude 
exhibits a well-pronounced minimum and resembles the 
theoretical dependencies in Fig|3] The effect is well pro- 
nounced though this local minimum corresponds to larger 
dispersion than in the noise-free case. In reality, the level 
of noise cannot be zero and this local minimum can pro- 
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Fig. 4: (Color online) Noise-induced regularization of chaotic 
dynamics in the Lorenz system at (a = 1, 6 = 8/3, r = 24.74) 
and (cr = 4, b = 1, r = 16): (a) variance of (z) calculated over 
20-loop segments of the chaotic orbit and (b) phase diffusion 
coefficient measuring coherence of chaotic oscillations. Noise 
amplitude a in graphs is scaled to the coordinate of the saddle 
points zq = r — 1. 



vide a smaller dispersion than that for the minimal noise 
level if, for instance, noise cannot be diminished below 
a = 0.05 zo (zq = r — 1 is the coordinate of the saddle 
points). Furthermore, the local minimum can become a 
global one for a larger inhomogeneity of the attractor. In 
Fig-HJa), the dependence of var(z) on the noise ampli- 
tude for a = 4, b = 1, r = 16 (where the chaotic attrac- 
tor nearly touches the slightly unstable non-trivial saddle 
points) has the minimum which provides 3 times smaller 
dispersion than the noise-free case. 

The dispersion of (z) calculated over finite segments of 
the chaotic trajectory is only a sample of a quantifier with 
which one can observe the ordering effect of noise we re- 
port. This ordering can be observed with other quantifiers 
of the system dynamics as well. For instance, one of the 
very important characteristics of the chaotic systems is the 
coherence of their oscillations. In particular, coherence de- 
termines susceptibility of the system to control forcing and 
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Fig. 5: (Color online) Noise-induced regularization of chaotic 
dynamics for mapping (0. 

predisposition to synchronization (e.g., see [19l[20]). The 
coherence is quantified by the diffusion coefficient of the 
chaotic oscillation phase (the oscillation phase to be intro- 
duced so that it growths by 2ir for one revolution of the 
trajectory). In Fig.@Jb) the dependencies of the phase dif- 
fusion coefficient on the noise amplitude are plotted for the 
same sets of parameters as we used for calculation of the 
dispersion of (z). One can see well pronounced minima. 
For a = 4, b = 1, r = 16 the diffusion can be suppressed 
by factor 5 in comparison with the noise-free case. 

Observation of the noise-induced regularization 
effect with mapping. — Our analytical treatment is 
not restricted to the Lorenz system which is utilized only 
for demonstration. The Lorenz system is known to possess 
significant peculiarities [18]; therefore it would be appro- 
priate to validate the effect we observe with independent 
examples. The classical scenario of the transition to chaos 
in the Lorenz system is known to be a special case of the 
transition to chaos via cascade of bifurcations of homo- 
clinic loops [5] . Within the vicinity of the bifurcation ac- 
cumulation point, a general system, where such a scenario 
of the transition to chaos occurs, can be described by the 
following mapping (in the noise- free case) [8]; 

J x v — pi , for x n > ; 
Xn+1 { -A 2 \x\ v +p2, for *„<(); 

, _ . -, / T\ ln(pi/a: n ) , for x n > ; 
tn+i - t„ + i + j T 2 \n(p 2 /\x n \), {oix n <0; 

where x n is the system state at the time moment t n , v is 
the saddle index, pi are the parameters featuring deviation 
from the bifurcation accumulation point (pi = p 2 = at 
this point). The iteration n — > n + 1 corresponds to the 
(n + l)-th revolution of the system — the oscillation phase 
grows for 2tt as n — > n + 1. For the sake of definiteness 
we restrict ourselves to the case of p\ = p 2 = p, v = 1/2, 



A 2 = \,T\ =T 2 = \ and introduce noise: 

x n+x = signum(x„)(v / Rd 7 - p) + a(n ; 

( 5 ) 

t n +i =t n + l + ln(p/\x n \). 

Here a is the noise amplitude and f„ are independent ran- 
dom numbers uniformly distributed in [— ^/3, y/3\ (i.e., 

(C 2 ) = l). 

For this mapping the phase grows uniformly with each 
iteration while time increments (t n+ i — t n ) are non- 
constant. The phase diffusion can be calculated from t n 
as a function of n (see, e.g., [H]). In Fig. [5] one can see 
that the regularizing effect of noise can be as well observed 
near p = 0.25 where the natural measure distribution is 
highly inhomogeneous. 

Conclusion. — We have reported the ordering effect 
of noise on the chaotic dynamics. This effect is rooted 
into inhomogeneity of the natural measure over the attrac- 
tor and the fact that noise diminishes this inhomogeneity. 
Without noise, the measure is determined by the largest 
multipliers of unstable periodic orbits (UPOs) composing 
the chaotic attractor and is typically extremely nonuni- 
form for long orbits due to the exponential dependence 
of the multiplier on the orbit length and the Lyapunov 
exponent which vary for different UPOs. Generally, the 
unweighted probability density function of some average 
for UPOs is centered around the value which does not have 
to correspond to the maximal measure (see Fig.[T] where 
the big tail of the weighted PDF clearly indicates that 
the natural measure of UPOs with large (z) is drastically 
larger than that of UPOs near the peak of the unweighted 
PDF). Hence, the dynamics is more contributed by UPOs 
which are not "geometrically typical" . Noise smooths the 
measure over the set of UPOs and reduces the role of ex- 
cursions along the "untypical" UPOs. On these grounds, 
we expect the effect to be a universal feature for chaotic 
dynamics with exception for some model systems with 
absolutely uniform natural measure (like the tent map). 
With strong enough inhomogeneity given, the ordering ef- 
fect can overwhelm the distorting action of noise on UPOs 
and lead to significant observable ordering of the system 
dynamics. We have developed the analytical theory for 
the effect of noise on the measure distribution over chaotic 
attractor and observed the ordering effects for dispersion 
of averages over finite segments of the chaotic trajectory 
and coherence of chaotic oscillations. Noteworthy, the re- 
ported phenomenon is novel; its mechanism being robust 
and vital is not related to any of well-known phenomena 
of noise-induced ordering (stochastic |T] or coherence res- 
onance (2] , suppression of deterministic chaos by noise [3] , 
etc.). 

* * * 
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